function I=gausshermi(f, n)
%I=gausshermi(f,n)
%用Gauss-Hermite 方法计算数值近似积分。
%

%Hermite 多项式
p=hermipol(n);
%多项式的根
x=roots(p(n+1,:));

   G=feval(f,x);		%节点上的函数值

%加权系数
for i=1:n
   C(i)=(2.^(n-1)*(factorial(n)).*sqrt(pi))./(n.^2.*(polyval(p(n,1:n),x(i))).^2);
end


I=dot(C,G);


end
